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Abstract 



Observations of binary inspiral in a single interferometric gravitational wave 
detector can be cataloged according to signal-to-noise ratio p and chirp mass 
^A. The distribution of events in a catalog composed of observations with p 
greater than a threshold po depends on the Hubble expansion, deceleration 
parameter, and cosmological constant, as well as the distribution of com- 
ponent masses in binary systems and evolutionary effects. In this paper I 
find general expressions, valid in any homogeneous and isotropic cosmological 
model, for the distribution with p and ^A of cataloged events; I also evalu- 
ate these distributions explicitly for relevant matter-dominated Priedmann- 
Robertson- Walker models and simple models of the neutron star mass dis- 
tribution. In matter dominated Friedmann-Robertson- Walker cosmological 
models advanced LIGO detectors will observe binary neutron star inspiral 
events with p > 8 from distances not exceeding approximately 2 Gpc, corre- 
sponding to redshifts of 0.48 (0.26) for h = 0.8 (0.5), at an estimated rate of 
1 per week. As the binary system mass increases so does the distance it can 
be seen, up to a limit: in a matter dominated Einstein-deSitter cosmological 
model with h = 0.8 (0.5) that limit is approximately z = 2.7 (1.7) for binaries 
consisting of two 10 Mq black holes. Cosmological tests based on catalogs of 
the kind discussed here depend on the distribution of cataloged events with 
p and A4. The distributions found here will play a pivotal role in testing 
cosmological models against our own universe and in constructing templates 
for the detection of cosmological inspiraling binary neutron stars and black 
holes. 

PACS numbers: 04.80.Nn, 04.30.Db 97.80.-d, 98.80.Es, 
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I. INTRODUCTION 



A. Overview 

The most promising anticipated source for the United States Laser Interferometer 
Gravitational-wave Observatory (LIGO) [Q], or its French/Itahan counterpart VIRGO is 
the radiation emitted during the final moments of inspiral before the coalescence of a neu- 
tron star - neutron star (ns-ns) binary system The instruments operating in both the 
LIGO and VIRGO facilities will evolve over time, eventually becoming sensitive to neutron 
star binary inspirals at distances approaching 2 Gpc 0]. 

Binary inspiral observations in the LIGO or VIRGO detectors will be characterized by 
their signal strength and "chirp mass" (a combination of the binary's component masses 
and cosmological redshift). The distribution of observed inspirals with signal strength and 
chirp mass depends on cosmological parameters that describe our universe (Hubble constant, 
deceleration parameter, density parameter), the distribution of neutron star masses in binary 
systems, the overall density of coalescing binaries, and the properties of the detector. In this 
paper I explore the binary inspiral event distribution (with signal strength and chirp mass) 
in the LIGO and VIRGO detectors for different cosmological models. 

In addition to their value as quantitative expectations of what LIGO and VIRGO can 
expect to observe, these a priori distributions will play a pivotal role both in the construction 
of templates for detecting binary inspirals and in the interpretation of the observations. The 
distributions presented here, evaluated for our preconceived notion of the binary inspiral 
rate and cosmological parameters, are the prior probabilities required to form the likelihood 
function from the observed detector output Q]. Additionally, the observed inspirals will be 
a sample drawn from a particular cosmological model characterized by Hq, qq, Qq, neutron 
star mass distribution, and evolution characteristic of our own universe. By comparing the 
observed distributions to the ones described here we can measure those properties of our 
own universe^] 

These cosmological tests are analogous to the number-count tests of classical cosmology, 
which, in their simplest form, involve observing the distribution of a source population 
as a function of apparent luminosity or redshift. The first suggestion that binary inspiral 
number counts be used to measure interesting cosmological parameters was made by Finn 
and Chernoff [||,|[ (although they did not use the language usually associated with this 
technique of classical astronomy). Using Monte Carlo simulations they demonstrated that 
the distribution of inspiral events with signal strength and chirp mass could be used to 
measure the Hubble constant. In this work I provide a more general, explicit, and complete 
exposition of the properties of binary inspiral observation catalogs. 



detailed study of how accurately those measurements can be made as a function of the number 
of inspiral observations is underway and will be published separately; here I focus on describing 
the properties of a catalog of observations defined by a data cut on the signal-to-noise ratio and 
how those properties can be used to measure cosmological parameters. 
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The cosmological implications of gravitational wave observations of binary inspiral were 
first recognized by Schutz 0. He pointed out that each inspiraling binary is a standard 
candle in the sense that, if observed in three independent interferometers, its luminosity 
distance can be determined from the observed detector response. If an observed inspiral is 
associated with one of the several galaxy clusters that reside in its positional error box (whose 
determination also requires three interferometers), and if the redshifts of those clusters 
are determined optically, then observation of several inspiraling binaries would lead to a 
statistical determination of the Hubble constant that is independent of the cosmic distance 
ladder and the uncertainties that lurk therein. 

Markovic proposed a variation on the general theme introduced by Schutz: he observed 
that known neutron star masses were all close to 1.4 and that, in any event, there is 
a maximum neutron star mass. The observed chirp mass is a function of the mass of the 
binary's two components and its redshift. Assuming that the mass distribution in neutron 
star binaries does not evolve significantly over the range of binary inspiral observations, 
examination of the chirp mass distribution in binary systems at fixed luminosity distance 
would reveal the corresponding redshift. Thus, gravitational radiation observations alone 
might suffice to determine the Hubble constant. 

Unfortunately, detailed calculations show that, even for the most advanced LIGO and 
VIRGO detectors that have been discussed, the fractional uncertainty in the measured lu- 
minosity distance will be of order unity for events seen more frequently than thrice per year 
{i.e., for events at distances greater than approximately 100 Mpc) and the angular po- 
sition error boxes for these events are likewise large (on order 10 deg^ [HOD- Consequently, 
cosmological tests that rely on accurate and precise measurements of the distance and po- 
sition of inspiraling binaries using LIGO and VIRGO are not promising. 

In contrast, the cosmological tests discussed here and in |^ require only gravitational 
wave observation in a single interferometer. Furthermore, advanced LIGO detectors can 
expect to observe approximately 50 ns-ns binary inspiral events per year, from distances up 
to 2 Gpc, whose signal strength can be measured to better than 10% and whose chirp mass 
can be measured to better than 0.1% The rate, depth, accuracy and precision of these 
single interferometer observations suggest that cosmological tests based on the distribution 
of observed events with signal strength and chirp mass have great promise. 



B. Outline 

In this paper I calculate the expected properties of a catalog of binary inspiral observa- 
tions made by a single interferometric gravitational wave detector. A catalog is presumed 
to contain a record of all binaries that coalesced during the observation period and whose 
inspiral signal-to- noise ratio p was greater than the catalog limit po- The catalog properties 
depend on the coalescence rate density and the binary system component mass distribution, 
both of which may vary with redshift. Notation for these and a discussion of the observations 
that bear on them is the topic of section 



In section |TI| I discuss the signal-to-noise ratio p as a measure of signal strength. The 
signal-to-noise ratio of a particular inspiraling binary depends on the binary's intrinsic prop- 
erties, its distance from and orientation with respect to the detector, and the detector's 
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intrinsic properties. The way in which p depends on the signal and detector properties 
suggests a useful measure of the detector bandwidth, which I discuss here as well. Finally, 
section concludes by specializing the discussion to the specific properties of the proposed 
LIGO and VIRGO gravitational radiation detectors. 

The important properties of a binary inspiral observation catalog are the distribution of 
cataloged events with p and the binary system "chirp mass" M.. These distributions depend 
on the cosmological model, which includes the evolution with redshift of the neutron star 
mass distribution and the coalescing binary number density. In section |^ I give general 
expression for these distributions and discuss how they may be used together with actual 
observations to test cosmological models. Also in this section I give expressions for the 
catalog depth (the maximum redshift of a binary system that can have p greater than po) 
and the total rate that proposed interferometers can expect to observe inspiraling binaries 
with p greater than po- All of these general expressions are evaluated explicitly for relevant 
matter-dominated Friedmann- Robertson- Walker cosmological models and a simple model of 
the neutron star mass distribution. Finally, I summarize my conclusions in section 0. 



II. COALESCENCE RATE DENSITY 



A. Introduction 



The signal-to-noise ratio p of a binary inspiral in a LIGO-like interferometer depends on 
the relationship between the binary and the detector [i.e., orientation, distance and redshift) 
and also on certain intrinsic properties of the system {i.e., component masses and spins). 
Of these intrinsic properties, the intrinsic chirp mass, 

Mo = p!"''M^"' (2.1) 

where p and M are the binary's reduced and total mass, plays the most important role: all 
the other intrinsic properties offer only small corrections to p. 

Gravitational-wave detectors like LIGO or VIRGO do not measure M.q] instead, they 
measure 



M = Mo{l + z), (2.2) 

where z is the system's redshift with respect to the detector. To distinguish between M., 
which involves the system's redshift, and M.o, which depends only on the binary's intrinsic 
properties, I refer to the former as the observed chirp mass, or simply the chirp mass, and 
the latter as the intrinsic chirp mass. 

In order to describe the binaries included in a signal-to-noise limited catalog we must 
first describe the coalescing binary distribution in space and in M.q. The notation I use to 
describe this distribution is defined in section |1I B| , while in section |II C| I discuss what is 
known about the distribution from present-day observations. 
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B. Definitions and notation 



Assume that coalescing binaries are distributed homogeneously and isotropically with 
the cosmological fluid and deflne the binary coalescence local specific rate density M by 

N= , (2.3) 

dtdVdMo 

where dV is a co-moving cosmological fluid volume element and dt is a proper time inter- 
val measured in the fluid rest frame. The total co-moving rate density on the surface of 
homogeneity at redshift z is thus 

n{z) = J X{Mo,z)dMo. (2.4) 

Deflne the ratio of the total co-moving rate density at a redshift z [n{z)] to that at the 
present epoch [no = n{z — 0)] by S: 

h{z) = S{z)ho. (2.5) 

The distribution of coalescing binaries with intrinsic chirp mass Jv[q on the surface at redshift 
z is 

P(>,„|,)=^^, (2.6) 

where, by construction, 

' dMoV{Mo\z). (2.7) 

The homogeneous and isotropic local speciflc rate density can thus be written 

A^(Mo, z) = h^S{z)V{M^\z). (2.8) 

Since we have deflned M and n in terms of the co-moving volume element and intrinsic 
chirp mass, in the absence of evolution £ is unity and V is independent of z. Additionally, 
for inflnitesimal dM.Q, V{M.Q\z)dM.Q can be interpreted as the probability that a randomly 
chosen binary on the surface at redshift z has intrinsic chirp mass in the range Alo to 
Mo + dM-Q. 

C. Observational constraints on H 

1. Rate density at the current epoch: ho 

Cosmological tests that depend on the observed distribution of inspirals with p and/or 
A4 do not depend on fiQ. Nevertheless, it is necessary to know fio in order to estimate how 
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long it will take to accumulate a catalog of observations large enough that such tests will 
give meaningful results. 

The best current estimate of the ns-ns binary coalescence rate density at the current 
epoch is 1.1 X 10~^hMpc~^ yr~^, where h is the Hubble constant measured in units of 
100 Kms~^ Mpc^^ [|ri|,|l^]. This estimates relies on the 3 observed binary pulsar systems 
that will coalesce in less than a Hubble time (PSRs 1913+16, 1534+12, and 2127+llC 
|T^-|T3|). Since the number of observed systems is small the actual rate is quite uncertain: 
Phinney ||12[ has estimated that, while unlikely, a rate two orders of magnitude higher or 
lower could be reconciled with current observations. 

Black hole - black hole (bh-bh) and bh-ns binaries are believed to form at rates compa- 
rable to ns-ns binaries; however, the masses of these systems and the fraction that merge in 



less than a Hubble time are entirely uncertain |]TT],|T2 



2. Intrinsic chirp mass distribution and evolution: V and E 

The intrinsic chirp mass distribution V{^M.q\z) depends on the binary system component 
mass distribution on the slice of homogeneity at redshift z. Denote a binary's component 
masses as mi and m2 and write their joint probability density on a surface of redshift z as 
P(mx^m2\z\ The intrinsic chirp mass distribution V{^M.q\z^ on that slice is then 



(2.9) 



The determination of V{^M.q\z^ thus reduces to finding P{m\^m2\z\ 

Both theoretical and observational evidence suggest that the neutron star mass distribu- 
tion is narrow [|T^-[TB|. A simple model of the mass distribution has the component masses 
in a binary uncorrelated and uniformly distributed between upper and lower bounds m„ and 
mi] then 



P{mi,m2\mi,mu) = P{mi\mi,my)P{ni2\mi,mu) 
= (m„ - nii)''^ 



(2.10) 



and 



P{Mo\mi,mu) 



rriu rm 



mi J mi 



drrii dm2- 



3 3 \ 1/5 



m\+m2 ) 



mi) 



(2.11) 



where mi and m„ may depend on z. The probability density P{^AQ\ml,mu) is maximized 
when A^o = {f^imuY^^ / {mi + m^)^/^. 

The limited observations of the neutron star mass distribution in binary pulsar systems 



provide independent 95% confidence intervals for m; and m„ [p!7|] : 
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1.01 < m^/Mo < 1.34 

1.43 < m„/M0 < 1.64. ^ ' ^ 



The most likely values of mi and are 

mi = 1.29 Mo, 

m„ = 1.45 Mq. ^^"^"^^ 

Over this narrow range P{Aio\mi,mu) is, to an excellent approximation, piecewise linear: 

P{Mo\mi,mu) ^ 

( MonH^ if mo > Mo> m<, 

{ if m> > A^o > mo, (2.14) 

I otherwise 



where 



m< = m//2^/^ (2.15a) 

"^0 = (2.15b) 

[mi + m„) 

m> = m„/2^/^ (2.15c) 

At present we can observe only nearby binary pulsar systems; consequently, there are 
no observations that bear directly on the variation oi £ oi V with z. Theoretical studies 
suggest that the initial mass of neutron stars formed by stellar core collapse do not vary 
significantly with the progenitor mass or composition [l^. After formation the mass may 
evolve owing to accretion from a companion; however, in any event it is not likely that 



either V or £ vary with z more rapidly than do galaxies. In section |^ I provide general 
expressions for and detailed examples of the expected distribution of events in a catalog of 
binary inspiral events; for the detailed examples I neglect evolution in £ and V entirely. As 
shown in section |IV D| , advanced LIGO and VIRGO detectors will observe ns-ns binaries 
from redshifts not expected to exceed z ~ 0.5 and with the preponderance of events arising 
from z ~ 0.1; consequently, neglect of evolution is not an unreasonable approximation. 



III. THE SIGNAL-TO-NOISE RATIO 
A. Outline 

The signal-to-noise ratio measures the signal amplitude in terms of a detector's noise 



properties. In subsection ^11 B| I define the signal-to-noise ratio p and discuss the subtle 
issue of how p is estimated, but not determined, by observation. The functional form of p 
depends on the detector response to binary inspiral radiation, which I describe in subsection 
p.11 C| . In subsection |111D| I give the binary inspiral signal-to-noise ratio in terms of the 
same parameters that characterize the detector response. The form of p suggests a natural 
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definition of a detector's effective bandwidtli for binary inspiral observations; I discuss tliis 
bandwidtli function in subsection [IIP as well. 

The relative orientation of the detector and binary is described by a function of four 
angles. While these angles cannot be measured by observations in a single interferometer, 
important properties of the angular orientation function can be calculated independently of 
the particulars of the binary or the detector. These properties play an important role in 
determining the binary inspiral catalog properties and interpreting individual observations. 



I discuss the angular orientation function in subsection HIE 



Finally, in subsection |TTTT] I specialize the discussion of the signal-to- noise to the proposed 
LIGO and VIRGO gravitational radiation detectors. 



B. The signal-to-noise ratio: general comments 

Let s{t) be the detector response to a gravitational radiation signal from any source. If 
the detector noise is Gaussian with one-sided power spectral density Sn{f) then the signal- 
to-noise ratio is defined to be 

where 's is the Fourier transform of the detector response, 

POO 

/ dfe'-^f'sit), (3.2) 



— oo 



and 's* is its complex conjugate |TU|. 

An observation of a gravitational wave signal in a noisy detector entails a measurement of 
the signal properties in the presence of a particular instance of the detector noise. Analysis 
of the detector output results in an estimate of p as well as other parameters that describe 
s{t). Throughout this work I use p to mean the actual signal-to-noise ratio, as defined by 
equation |5'. i|, and not the estimate that arises in an observation. 



In section |I^ I find the distribution of sources with p and Ai in different cosmological 
models. By comparing these distributions with the observed one we can determine the model 
that best describes our own universe. In making that comparison it is critical to distinguish 



between p and (as defined by equations |3.1| and |2.2D and the estimates of p and At that 



results from observations made in a detector. 

The estimate that results from an observation is a probability distribution P for the 
parameters that describe the signal — in the case of binary inspiral in a single interferometer, 
these include p and the chirp mass Ai (see sec. [Ill Q below). The probability distribution 
associated with an observation is generally not reported; instead, what is reported most 
often is a set of estimators that characterize the distribution and its moments. Among the 
most popular is the the maximum likelihood estimator, which is the set of parameter values 
that maximize P. 

Estimators are summaries of the distribution P and their utility depends on how ac- 
curately they are able to represent it. When P is very sharply peaked {i.e., there is little 
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uncertainty in the measurement) then it may be approximated by a 5-function and summa- 
rized accurately by the maximum hkehhood estimator. When P is sharply peaked but with 
a not insignificant width, then it may be approximated near its peak by a Gaussian and ac- 
curately represented by the maximum likelihood estimators and their covariance. When the 
distribution is not sharply peaked, however, then the more general uncertainties refiected in 
the detailed structure of P play an essential role in the observation's interpretation and no 
summary of P is especially useful. 

That a small set of estimators built from P for a particular observation do not provide 
a useful summary does not mean that the observation itself is unreliable or uninformative; 



rather, it means only that greater care must be taken in its interpretation. Finn ||2y, sec. 
(c)] gives an example of how the the probability distribution resulting from an observation 
should be used in the interpretation of neutron star mass observations in binary pulsar 
systems; a further discussion of this point in the context of cosmological tests using binary 
inspiral observation catalogs is part of a work in preparation. 



C. Detector response to binary inspiral 

The detector response s{t) to the gravitational radiation from an inspiraling binary sys- 
tem depends on the distance and relative orientation of the source and the detector, as well 
as on certain intrinsic properties of the binary. The relative orientation of the source and 
the detector is described by four angles: two {9 and 0) describe the direction to the binary 
relative to the detector, and two (i and ip) describe the binary's orientation relative to the 
line-of-sight between it and the detector. 

To describe 6 and 0, consider a single interferometric gravitational wave detector whose 
arms form a right angle. Let the arms themselves determine the x and y axes of a right- 
handed Cartesian coordinate system with the z-axis pointing skyward. In this coordinate 
system the gravitational waves from an inspiraling binary arrive from the direction n, which 
can be defined in terms of the spherical coordinates 6 and in the usual way: 

cos 6* = —n ■ z, (3.3a) 

tan0= —. (3.3b) 

n ■ X 

To describe i and ip, let J represent the total angular momentum of the binary system. 
The detector response s{t) depends on the inclination angle i between J and n, 

cosi = —J ■ n/\J\, (3.3c) 

and the orientation ip of the angular momentum about n, 

cot^ = ^ ^ ^ ^ — -T. (3.3d) 
J -yz — n[z ■ n)\ 

The conventions for the orientation angles 6, (p, i, and ip described here are the same as 
those used in MM (note, however, that the description of the angles in M is incorrect). 
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At Newtonian [i.e., quadrupole formula) order the only intrinsic property of the binary 
system that affects the waveform is Aio. At this order the detector response (a dimensionless 
strain) to the binary inspiral signal is 

{irfMf^^ cos [x + $(t)] for t < T', 

for t > T\ 



s{t) 



^0 




(3.4a) 



where x is a constant phase, 



e = 



+ AF^ COS- I 



Fl (l + cos2i) 
(l + cos^ 6) cos 2(j) cos 2ip 



11/2 



/ 



TlM 



— COS 6 sin 20 sin 2-0, 
-f cos^ 9) cos 20 sin 2tp 
+ cos 9 sin 20 cos 2'0, 
5 M 



256T-t 



3/8 



<^ = 2n f{t)dt = -2 



T 



T-t 
5M 



5/8 



(3.4b) 

(3.4c) 

(3.4d) 
(3.4e) 

(3.4f) 



di is the source's luminosity distance, T' is the moment when the inspiral waveform termi- 
nates (either because the binary components have coalesced or because the orbital evolution 
is no longer adiabatic), and T > T' would be the moment of coalescence if the two compo- 
nents of the binary system were treated as point masses in the quadrupole approximation 
(the difference between T and T' is small but not negligible). Note that dependence of the 
response on the angles 9, 0, i, and is contained entirely in the orientation function 0, 
which does not depend on any other properties of the binary or the interferometric detector. 
This important point will be discussed further in subsection |III E . 

Post-Newtonian corrections do not contribute significantly to the signal-to-noise ratio for 
solar-mass binary inspiral in the LIGO and VIRGO interferometers. For symmetric binaries 
[i.e., those with equal or near equal mass components) the first post-Newtonian correction 
is proportional to M/r, where M is the binary's total mass and r the component separation. 
Advanced LIGO and VIRGO interferometers are expected to be most sensitive to binary 
inspiral radiation in the bandwidth 20-200 Hz (see sec. p.11 D 2 ) and, in this bandwidth, 
M/r < 4% for solar mass binaries. Post-Newtonian effects are more important for more 
massive binaries, since at fixed quadrupole radiation frequency M/r is greater for greater 
M: 



r J 



0.042 
0.16 



/ M 
200 Hz2.8M0_ 

/ M ^'^'^ 



2/3 



200 Hz 20 Mq 



(3.5a) 
(3.5b) 
(3.5c) 
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Post-Newtonian effects are more important for binaries with extreme mass ratios: in general, 
the first post-Newtonian correction is proportional to {6m/ M){M /ry^^ , where Sm is the 
component mass difference. Finally, spin-orbit coupling in systems whose components have 
large spin angular momenta can lead to orbital precession and a waveform modulation which 
can affect p significantly ||21|| . 

That the quadrupole waveform can be used to estimate p for an inspiraling binary does 
not mean it can also be used to detect a binary inspiral signal in the detector output. 
Detection, or estimation of signal parameters from observation, involves the comparison 
of the detector output with a model of the detector response to the radiation. "Small" 
differences between the actual and modeled signal evolution, particularly its phase, can have 
a large effect on the measured value of the signal parameters, including the signal-to-noise 
ratio estimated from the observations P,|22[ . 



D. Binary inspiral signal-to- noise ratio 



The signal-to- noise ratio p (eq. ^.Ij ) corresponding to the detector response s{t) (eq. 



Oa]) is 



p = 8eli 



M 



1.2 



5/6 



(3.6a) 



where 



Tn = 



X7/3 
/max 



C(/n 



19271 V2oy 



5/3 



— 1 X7/3M2 



0' 



df {tiMqj 



(7r/Mo)'/^^„(/)' 
/(T')/2, 

1 /-^-^-^ df {ttUq^ 



X7/3 Jo 



(3.6b) 

(3.6c) 
(3.6d) 
(3.6e) 



In equation |3.6a| 0, A4 and di depend on the particular binary system under consideration 
and vq is a characteristic distance that depends only on the detector's noise power spectral 
density Sn{f)- The dimensionless function (, which also depends only on the detector's 
noise spectrum, increases monotonically from to 1. Its argument /max is the redshifted 
instantaneous orbital frequency when the inspiral terminates (at t = T') either because the 
compact components have coalesced or because the orbital evolution is no longer adiabatic 
and coalescence is imminent. 



1. ro and C, 



The characteristic distance tq and the function ( describe different aspects of an inter- 
ferometric detectors sensitivity to binary inspiral gravitational radiation. For a fixed binary. 
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the larger ro, the greater p; for fixed p, the larger ro, the farther the detector can "see." 
Decreasing the noise power in any band increases tq, but (owing to the factor /"^^^ in the 
integrand of the expression for X7/3) improvements at low frequency are more effective than 
those high frequency. 

The function ( reflects the overlap of the signal power with the detector bandwidth. The 
orbit of an inspiraling compact binary evolves adiabatically owing to gravitational radiation 
emission until an innermost circular orbit (ICO) is reached at an instantaneous orbital 
frequency /ico- At the ICO the orbit evolves on a dynamic timescale and the components 
coalesce quickly. Thus, the adiabatic inspiral waveform ends when the orbital frequency 
reaches /ico0- The quadrupole radiation frequency observed at the detector when the orbital 
frequency is /ico is 

2/max = ~r~ (3-7) 
-L i ^ 

and the inspiral signal termination is represented in the detector response by the cut-off at 
t > T' (see eqs. p.4a] and |3.6d|) . 



Since the binary orbit evolves adiabatically from low frequencies to /ico? the observed 
quadrupole radiation spectrum has signiflcant power only up to frequency 2/max- If 2/max is 
very much greater than the frequency where the detector noise is minimized, than the signal 
power bandwidth overlaps completely with the detector bandwidth and C(/max) — 1- On the 
other hand, if 2/max is much less than the frequency where the detector noise is minimized, 
than the overlap of the signal power and detector bandwidths is negligible, C(/max) — and 
p ^ 0. 

The orbital frequency at the transition of the binary orbit from adiabatic inspiral to 
plunge and coalescence has been studied using high-order post-Newtonian methods p3[. For 



symmetric binaries {i.e., those with equal-mass components) the instantaneous redshifted 
instantaneous orbital frequency of the ICO is given by 

/max = z — ; — ) (3.8a) 



/ico 








710. Hz 


/2.8M 


l + z 


V M 


99. Hz / 


'20. M0 


l + z \ 


M 



(3.8b) 
(3.8c) 

where M is the binary's total mass. More generally, for binaries with compact components 
/ico depends inversely on M and directly on a function of the dimensionless ratio fJ^/M; 
Kidder ||2^, flg. 4] shows M/jco for binary systems of arbitrary mass asymmetry. For a 



^The radiation waveform from the final plunge and the early stages of coalescence is not yet 
known. If, after coalescence, a black hole forms, then the final radiation reflects the black hole's 
quasinormal modes, which are damped rapidly. The inspiral waveform, and with it our ability to 
model the detector response, ends when the orbital frequency reaches /ico- 
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ns-ns binary the component's proper separation at the ICO is greater than a neutron star 
diameter; so, coalescence occurs after the transition from inspiral to plunge [23|. Tidal 



dissipation is important in determining the inspiral rate only in the last few orbits before 
contact [^,^]; consequently, equation ^]8| is applicable for ns-ns binaries and should be 
applicable for black hole binaries as well. 



2. Detector bandwidth and data analysis templates 

The probability that a signal with detector response s{t) is present in the output g{t) of 
a noisy detector is related to 

2{g,s)-{s,s), (3.9a) 

where 

The actual signal, and consequently the detector response s, may be difficult or impossible 
to evaluate exactly (this is certainly the case for binary inspiral); so, we would like to be 
able to use an approximate model S{t) in lieu of the actual detector response s in equations 



From equations ^]9| it is clear that any approximate model S{t) may be used as long as 
it matches s'(/) closely wherever |sp/5'n(/) is relatively large. This latter quantity is just 
the fraction of contributed by the signal power at frequency / and is proportional to 
This suggests that we define the bandwidth function 

Bif) = { 2C(/,na.) ^ ^ (3.10) 



which is the fraction of the signal-to-noise ratio contributed by signal power at frequency / 
in a logarithmic frequency bandwidth. An approximate detector response model S may be 
used instead of the actual response 's as long as S accurately reflects 's wherever B is large. 
Knowing where and how S needs to be accurate can simplify greatly the construction of 
approximate templates for data analysis; for this purpose B should prove a useful guide. 



E. Distribution of the orientation function Q 

The signal-to-noise ratio of an inspiraling binary in a single interferometric detector 
depends on the relative orientation of the source and the detector through the angles 9, 0, i, 
and "0. The dependence of p on these angles is confined to the angular orientation function 
O, which is independent of all other properties of the binary system. From observations of 
binary inspiral in a single interferometer we can measure p and A4 but not G. Even though 
9 cannot be measured, because it depends only on 9, 0, i and ip we actually know quite a bit 
about it: since, averaged over many binaries, cos 6', (p/n, cosi, and ip/ir are all uncorrelated 
and uniformly distributed over the range (—1, 1) we know the probability that takes on 
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any value. This probability distribution is found numerically in 
< < 4 and that to an excellent approximation 



Pe(e) 



56 (4 - 6)^ /256 if < < 4 
otherwise. 



here I note only that 



(3.11) 



To determine the binary coalescence rate p greater than po we need to know M and also 
how the signal-to-noise ratio of binaries with intrinsic chirp mass M.QO\ia. surface at redshift 
z is distributed. This latter distribution, Pp{p\M.Q^ z), is related to Pq: 



PMMo.z,C,V) = Pe[e{p)] 



Pr 



dQ 
dp 



Mo,z 



1.2 







X 



roil + zf V M 



dL{1.2MQ/Mo 



5/6' 



,5/6 



^ro (1 + Z 



5/6 



(3.12) 



where C represents the cosmological model (Hq, qq, Qq, V, £) and V the detector [i.e., tq 
and Q. 



F. ro and C for the LIGO and VIRGO detectors 

LIGO will consist of three interferometers: one in Livingston, Louisiana and two in 
Hanford, Washington. The Louisiana interferometer and one of the Washington interferom- 
eters will have 4 Km armlengths; the second Washington interferometer will have a 2 Km 
armlength and share the same vacuum system as the 4 Km interferometer. For the pro- 
posed LIGO interferometers as described in and modeled in , tq ranges from 13 Mpc 



(for "initial" interferometers) to 237 Mpc (for "advanced" interferometers)^. As the LIGO 
detectors develop, incremental improvements will increase tq. 

The orientation of the Washington and Louisiana interferometers were chosen to be as 
close to parallel as possible; consequently, a simple approximation treats the LIGO detector 
(all three interferometers operating in "triple-coincidence") as a single interferometer of arm 
length 

[42 + 42 + 2^] ^'"^ Km = 6 Km. (3.13) 
For this "super-interferometer" tq ranges from 19 Mpc to 355 Mpc. 



^The noise model used in [|2^ and assumed a thermal noise spectrum corresponding to vis- 
cous damping in the pendulum suspensions and the internal modes of the test masses. It is now 



recognized that these modes are structure damped [27-p9| with a correspondingly different noise 



spectrum. Preliminary estimates indicate that this improved noise model reduces tq for initial 
LIGO interferometers, but leaves it unchanged for advanced ones. 
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Here and below, reference to the LIGO detector refers to the three interferometers op- 
erating as a single detector, while reference to a LIGO interferometer indicates one of the 
4 Km interferometers operating in isolation. 

Figure |l] shows ( for early LIGO and VIRGO interferometers and also more advanced 
LIGO interferometers as have been discussed in the literature [^,0. The solid curve shows 
( representative of advanced interferometers, while the dashed and dotted curves show ^ 
representative of initial LIGO and VIRGO interferometers. As the LIGO detector evolves, rg 
will increase from approximately 20 Mpc toward 350 Mpc, and ( will evolve from the dashed 
curve toward the solid curve. Setting aside the overall sensitivity gain with increasing tq, 
the evolution of ( signifies an increasing relative sensitivity to systems with small /max- ^-c, 
systems with larger total masses and/or redshifts. 

VIRGO will consist of a single 3 Km interferometer. The target noise curve for the 
initial operation of the VIRGO interferometer is described in ||^; for this interferometer 
tq is 13 Mpc. Like the LIGO interferometers, incremental improvements in VIRGO will 
increase tq and it is reasonable to assume that these improvements could raise tq for VIRGO 
to approximately 200 Mpc. The dotted curve in figure |1] shows C for the initial VIRGO 
interferometers. Note how, even though tq is the same for both the initial LIGO and 
VIRGO interferometers, the relative sensitivity of these two detectors to signals at low and 
high frequency — corresponding to more or less massive binaries — is markedly different. 
Improvements in the VIRGO interferometer can be expected to evolve ( toward the solid 
curve as well. Note that the initial VIRGO interferometers are expected to be relatively 
more sensitive to massive {i.e., low /max) binaries than are the initial LIGO interferometers. 

For neutron star binaries /max — 710 Hz (eq. |378| ) and C(/max) — 1 in any of the proposed 
LIGO and VIRGO interferometers (see fig. 0). Figure ^ shows the bandwidth function for 
neutron star binary observations in early VIRGO (dotted curve) and LIGO (dashed curve) 
interferometers and also more advanced LIGO interferometers (solid curve). It is clear that 
the approximate detector response models for binary inspiral observations must evolve with 
the interferometers: in the early interferometers the detector response models will need to 
be quite accurate at high frequencies (/ ~ 200 Hz) but not at low ones; on the other hand, 
in more advanced interferometers these models will need to be accurate at low frequencies 
(/ ^ 60 Hz) and, to obtain this low frequency accuracy, the high frequency performance 
may be sacrificed at no cost to the signal detectability. 



IV. DISCUSSION 
A. Introduction 

Consider a catalog of binary inspiral events with signal-to-noise ratio p greater than a 
threshold po- Suppose that for each event in the catalog p and Ai are known. How do we 
make use of the catalog to test cosmological models against cataloged observations? 

The simplest test involves the distribution of events with p. As the catalog limit po 
decreases, sources at increasingly larger distance become members. The number of added 
sources depends on the increase of the spatial volume, the density of sources at that distance, 
and (since these sources are events that occur at a given rate) the cosmological redshift. 
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Thus, adopting a cosmological model C implies an expected distribution P{p\pQ,C,V) for 
the catalog events taken with detector V. Denote the cataloged binary inspiral signal-to- 
noise ratio observations by {p\p > po}. Suppose that, before studying these observations, 
we have reason or are otherwise prejudiced to believe the probability that C is the correct 
cosmological model is P{C). Using Bayes' law of conditional probabilities, the a posteriori 
probability that we assign to model C after considering the observations is 

P{C\{p\p > Po}) oc P{C)l[P{p,\po,C,V,), (4.1) 

i 

where pi is the signal-to-noise ratio of the i^^ catalog observation and Vi represents the 
detector configuration (ro, C) when the observation was madeQ. This test is exactly analogous 
to number-flux cosmological tests using distant galaxies. 

More subtle tests involve other distributions of cataloged events. For example, each 
source in the catalog is characterized by its observed chirp mass Ai. The observed chirp mass 
depends on both the intrinsic chirp mass Aio and the redshift z (see eq. p.2| ); consequently, 
the distribution of events in the catalog with Ai depends on the cosmological model. Denote 
the cataloged chirp mass observations by {M.\p > po}. Adopting a cosmological model C 
implies an expected distribution P{Ai\po, C, T>). As before, if we initially favor model C with 
probability -P(C), then after studying the observations the probability we ascribe to model 
C is 

P{C\{M\p > Po}) oc P{C)l[P{Mi\po,C,Vi), (4.2) 

i 

where Aii represents the i^^ cataloged chirp mass observation. 

Of the four parameters that describe, at quadrupole order, binary inspiral observed 
in a single interferometer, only p and Ai convey astrophysically interesting information. 
The distributions P{p\pQ,C,V) and P{A4\pq,C,V) are each integrals over the distribution 
P{p, J^\po,C,V) that completely characterizes the catalog: 



P{p\po,C,V) = J P{p,M\po,C,V)dM (4.3a) 

/•oo 

P{M\po,C,V)= / P{p,M\po,C,V)dp. (4.3b) 

J po 

These integrals are summaries of the catalog contents: as such, they are less informative 
than P{p, Ai\po,C,V). The most sensitive test that we can make using the catalog involves 
not a summary, but the full information available: given the observations {p,Ai\p > po}, 
the probability that cosmological model C describes our universe is 



^Changes in the detector noise spectrum change the relation between p and the cosmological 
model. This is not a problem in the interpretation of a catalog as long as the detector properties 
are properly associated with each cataloged observation. 
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P(C|{p, M}, po) oc P(C) H P{p„ Mi\po, C, Vi). 

i 



(4.4) 



Cosmological tests based on the summary distributions are still useful to make. Summary 
distribution often depend only weakly or not at all on some of the parameters of the model 
C] in that case, the effective dimensionality of C is reduced in the summary test and we may 
be able to distinguish more closely among cosmological models described by the remaining 
parameters than with a test using the full distribution. This is particularly true when the 
number of observations in the catalog is small. 



B. Outline 

In this section I give general expressions for the distributions P{p, JiA\po,C,T>) and the 
summary distributions P(p|po; C, V) and P(A^|po, C, T>) that can be evaluated in any homo- 
geneous and isotropic cosmological model; in addition, I evaluate the summary distributions 
explicitly for matter-dominated Friedmann- Robertson- Walker (FRW) models. The essential 
properties of FRW cosmological models are summarized in section [IV (J| . 

In addition to the distributions P{p, A4\po,C,V), P{p\pQ,C,V) and P{J^\po,C,V) are 
other catalog properties of intrinsic interest. In a given cosmological model there is a distance 
beyond which no binary inspiral of fixed intrinsic chirp mass will have signal-to-noise ratio 
greater than the catalog limit po, and I calculate this catalog depth Zq in section P^VD] . 
A catalog takes time to build, and, in a given period of time, the size of a binary inspiral 
catalog is limited by the rate dN/ dt at which binaries coalesce with p greater than the catalog 
threshold po. In section [TV E| I calculate both the expected rate of binary inspiral observations 
and the distribution P(p|po, C, V) for advanced LIGO interferometers and describe how these 
depend on the cosmological model C. Finally, in section |1V b] 1 calculate the distribution 
P{J^\po,C,V) and describe how it depends on C. 



C. Cosmological model 



Specific examples of the catalog distributions, catalog depth, and rate of binary inspi- 
ral observations made in this section are in the context of matter-dominated Friedmann- 
Robertson- Walker (FRW) cosmological models. The Robertson- Walker spacetime metric 
has the line element 



ds' 



-dt^ + a^{t) 



dr^ 



+ {de^ + ^in^ 



(4.5) 



where a{t) is the usual (dimensioned) scale factor, r is a dimensionless parameter related 
to the area of spheres of constant radius, and k is —1, or depending on whether the 
spatial geometry of the slices of homogeneity has positive, negative, or zero curvature {i.e., is 
closed, open, or flat). In matter dominated FRW models (which have vanishing cosmological 
constant) the co-moving radial coordinate r and the luminosity distance di can be written 
in terms of the redshift explicitly: 
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zqQ + (go - 1) (1 + 2go2) 
r = 



1/2 



dL = aoil + z)r, (4.6b) 

where Qq, Hq, and go are the scale factor, expansion rate, and curvature parameter (decel- 
eration parameter) at the present epoch eq. 15.3.24]. 



D. Sample depth 

The signal-to-noise ratio of an inspiraling neutron star binary system with intrinsic chirp 
mass A^o is 

Since is between and 4, even an optimally oriented binary system has p less than po 
when z is greater than zq, where Zq satisfies 



PodLizo) /1.2M^^'/^ 



^roC(/max)^/'^ V Mo 



(4.8) 



Evaluation of zq requires C(/max), which depends on the details of the detector's noise power 
spectral density (through () as well as the binary's component masses and redshift (through 
/max) ■ For advanced interferometers C — 1 as long as /max ^ 70 Hz (see fig. |ID while /max — 
710 Hz [2.8 Mq/(1 -|- 2;o)M] for symmetric binary inspiral (see eq. |3.8| ). Consequently, we 
can approximate C — 1 as long as 1 -|- ^ 10 (2.8 Mq/M). For small zq we can approximate 
— z/Hq] then 

32goro f Mo 



Po V 1-2 Me 

^OASh(^)(^^)(^Y\ (4.9b) 
\poJ V355Mpc; V1-2M0; ^ ^ 

which is much less than 10. Thus, the approximation ~ 1 is a good one for binary systems 
with solar mass components, but not for binaries whose components are on order SM©. 
For a specific example, focus attention first on solar mass component binaries. Then 
~ 1 for any of the proposed LIGO or VIRGO interferometers (at the end of this subsection 
I return to consider briefly the case of more massive binary systems where this is not true). 
In an Einstein-deSitter (go = 1/2) cosmological model zq is then given explicitly by 



where 



■ 256 2 16 

/5+^«o + y«o 



1, (4.10a) 
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SHoro ( Mo 



Po 



1.2 M„ 



0.12/1 ( — 



5/6 



Mi 



5/6 



355Mpcy V1-2M0^ 
h is the Hubble parameter [Hq in units of 100 Km/s/Mpc), and 



(4.10b) 
(4.10c) 



P 



/leaoV , , 4096 3 



2 ■ V 3 ; +U + ^"° 



1/2- 



1/3 



More generally (go 7^ 1/2 but ( still unity) Zq satisfies 

„6 1 



2:0 = X 



(4.10d) 



(4.11a) 



where a; is a root of 



= x'^'^ + SqoaoX^^ + IGqQaQX 

- 2qoX^ - 8ao(2go - + 2qo-l. 

The appropriate root of this equation can be found as a power series for small ao- 
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(4.11b) 



4«o + - (3go + 2) ^ 8 (6go - 1) 



+ 



160 
81 



(54g2 + 9go - 7) + 



5600 
243 



27ql - ^qo + |^~) a'o + 0{al 



(4.12) 



This truncated expansion is accurate to better than 0.2% for ao < 0.12 and < go ^ 1; for 
ao < 0.1 and < go < 3/4 the accuracy is better than 0.01%. An asymptotic expansion for 
zq valid for large ao (but ( still 1) is 



4096g^a|j - ^uJ2ql (go - 1) a^ + 6go - 1 + O (ao"^) 



(4.13) 



This asymptotic expansion is accurate to better than 1% for ao > 1/2 and 0.1 < go ^ 1- 
Figure |^ shows as a function of h for neutron star binaries with M.q = 1.19 
(corresponding to 1.37 Mq neutron stars), po = 8; ''^0 = 355, and three different values of go 
corresponding to an open (go = 1/4, dotted curve), flat (go = 1/2, sohd curve) and closed 
(go = 3/4, dashed curve) cosmological model. The redshift of the most distant such source 
is less than 1/2 for h < 0.8 and does not exceed 7/10 for h < 1. The sensitivity of Zq to go 
is modest but significant: inspection of equation [4.11| shows that 



4ao [1 + ao (2go - 1)] + O [a^ (2go - 1)'] 



(4.14) 



where ao — 0.12h for advanced LIGO interferometers. For open spatial geometries (2go — 1 < 
0) the most distant sources are at smaller redshifts than in closed spatial geometries (where 
2go - 1 > 0). 

In luminosity distance, the sample depth is 
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_ 4ao 



40 . 40 ... 2 , 1280 ,^ 3 



1 + y «o + y (2go + 1) c^o + ^ (9^0 - 2) «^ 

(4.15) 



40 

+ ^ (2052g2 - i80go - 143) + O (a*]) 



Figure @ shows (iL,o for the same cases {h, go) as figure || shows zq. Advanced LIGO inter- 
ferometers may observe neutron star binaries with p greater than 8 at luminosity distances 
of order 2 Gpc. 

For more massive binaries C(/max) is substantially less than unity at zq and the approx- 
imation ^ 1 is no longer valid. Figure |^ shows zq (given implicitly by eq. |4.8| ) for an 
advanced LIGO detector (ro = 355 Mpc) and a range of binary systems in two distinct 
Einstein-de Sitter (go = 1/2) cosmological models. One pair of curves show zq for symmet- 
ric binaries consisting of two {e.g.) black holes each of mass Mbh, while the second pair 
of curves shows for asymmetric binaries, consisting of {e.g.,) a black hole of mass Mbh 
and a neutron star of mass 1.4 Mq. Note how, in all cases, zq does not increase without 
bound as Mbh increases; instead, it has a maximum value for some Mbh and decreases as 
Mbh increases further. 

The finite bandwidth of the detector and its overlap with the signal power, as represented 
by C, determines the maximum redshift at which a detector can observe binary inspiral. If 
^ were constant then zq would increase monotonically with Jv[q (see eq. ^Tsj) ; however, C, is 
not constant and as M. increases /max decreases both because /ico is decreasing and also 
because Zq is increasing (/max oc M/(l + z); see eq. |3l8| ). Eventually C begins to decrease 
with increasing M.q or decreasing /max (for advanced LIGO interferometers, C, remains unity 
until /max — 70 Hz; see fig. [l|). Once C begins to decrease so to will zq. For sufficiently large 
M.Q even a local binary {z = 0) will coalesce before the gravitational radiation signal enters 
the detector bandwidth. 

Thus, decreasing the detector noise at low frequencies has a dual effect: by increasing 
vq it increases the overall detector sensitivity, and by decreasing the frequency below which 
C -C 1 it further boosts the sensitivity of the interferometer to more massive and more 
distant binary systems. Advanced LIGO detectors are expected to be most sensitive to 
inspiraling binaries with two 10 Mq black hole components, which will be observable with p 
greater than 8 at redshifts greater than 1.5. 



E. Coalescence rate above threshold 

The distribution of catalog events with p is given by 



where 



(fNjdtdp 

P{p\p,,C,V) = { /^M P>P^^ (4.16a) 
p < Po, 



dM,—r-^, (4.16b) 



dtdp J dt dp dj\4.{^ 
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dt dp dAii 



dz 



dr AiraQr"^ hQS{z) 



X 



dN 
'dt 



ViMo\z) 
1 + z 
d^N 



P,{p\z,Mo,C,V), 



dtdp' 



(4.16c) 
(4.16d) 



and Pp is given in equation |3.12 



Since d N/dtdp depends on p only through Pp, the integral over p in equation [4.16d| can 



be evaluated explicitly. First note that 

dpPp{p\Mo,z,C,V) 



PO 



Ce{x) 



where 



X 



Po 



dr 



8 ro (1 + z) 
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(4.17a) 
(4.17b) 

(4.17c) 



The probability density Pp is a conditional one: it depends on Aio, z, and the cosmological 
model. In contrast, the distribution Pq of G is universal: it is independent of all the specific 
properties of the binary or the detector. Cq{x) is universal in this same way and can be 
evaluated with Pq. In terms of Cq the total rate of inspirals with p greater than po is 



dN 
Itt 



— = dM 



d'^N 
^dtdMo'' 



(4.18a) 



where 



d^N 
dtdMo 







dz vT 



kr'^ I + z 
dAz) 



ttO (1 + ^)^^^C(/max) 



(4.18b) 



and 



8Foro / Mq 



Po 



1.2 Mq 
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(4.18c) 



Equation |4.18| for dN/dt, which is also the total rate of binary inspiral with p > po, and 
equations [4.16| for P{p\pQ,C,V) are valid for any homogeneous and isotropic cosmological 
model. 
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1. Examples: Constant Mq 



Turn now to the specific example of neutron star binaries in matter-dominated FRW 
cosmological models, so that ~ 1 and r and di are given by equations |4.6| . Assume also 
that evolution in the binary population is unimportant over redshifts z < 1/2; then S is 
unity and V is independent of z. Finally, use the approximation given above for Pq (eq. p. 11 
to evaluate Cq: 



Ce{x) 



1 if a; < 

(1 + a;)(4 - x)^/25Q if < x < 4 
if 4 < X. 



(4.19) 



The integral for (pN/dtdM.Q (eq. ^.18b| ) can then be evaluated as a power series in small 



d^N fdN\ , , , 



dtdMo \ dt J Q 



(4.20a) 



where 
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1 - —ao + — (6go + 85) al 
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-^qoal + ^ (I188g2 + 9108go - 4123) a', + O (a^) 



and cto is given by equation |4.10b| . Similarly, d^N/dtdpdJ^o is given by 



d^N 



dtdpdM.Q \dtdp J Q 



d^N 



V{Mo) 



(4.20b) 



(4.21a) 



where 
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dt dp 
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33 ^ 24057 
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(4.21b) 



and 
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(4.21c) 



Equations [4.20| and [4.21| for d^N/dt dMo and d^N/dt dp dMo can be integrated over A^o to 
find dN/dt and d'^N/dt dp, and thus P{p\po,C,V). 
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Observational and theoretical evidence suggests that the ns-ns binary intrinsic chirp 
mass distribution is sharply peaked (see sec. |II Q ). In the most extreme limit that all binary 
systems have intrinsic chirp mass A^o^ ^ is a (5-f unction in Mq and 



dN fdN\ _ , 

1 — )^ (4-22a) 



dt V d'l^ 



P{p\po.C,V) 



{d^N/dtdp)o 
{dN/dt)o 
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20 /'11QQ^2 I nins^_ _/110Q^^,4 



+ 0{ala^). (4.22b) 



Figure ^ shows {dN/dt)^ as a function of /i in an Einstein-deSitter cosmological model 
(go = 1/2) for a reasonable signal-to- noise ratio limit in an advanced LIGO detector (po = 8 
and ro = 355 Mpc), typical neutron star masses (A^o = 1.19 M©), and two different coales- 
cence rate densities. The solid curve shows {dN/dt)^ for the "best guess" coalescence rate 
density, no = 1.1 x 10~^/i Mpc~^ s~^ (see sec. [11 (J 1| ), which is proportional to h\ the dashed 
curve shows {dN/dt)o for a constant coalescence rate density no = 8 x 10"® Mpc~^ yr^^- The 
dashed curve shows how {dN/dt)Q scales for constant no in advanced interferometers with 
interesting h. Present estimates of no suggest that, ii h = 0.75, advanced LIGO detectors 
can expect to observe a little more than 1 neutron star binary inspiral event per week. 

The distribution P{p\pQ,C,V) is also sensitive to h. Let Co represent a flat and static 
cosmological model. Then 

Pip\Po,Co,V) = !^f^/p' P^^/^l (4.23) 

The first order correction (in 1/p) owing to the expansion of the universe depends only on 
the rate of expansion (Hq); corrections owing to the curvature of space {k ^ 0) enter only 
at second order in 1/p. Since, for ns-ns binaries and reasonable limiting signal-to-noise is 
less than 1% even for advanced interferometer designs, cosmological tests that focus only on 
the observed distribution of p will be insensitive to go even though the redshift to the most 
distant sources in the catalog is large (this weak dependence on go is characteristic of number- 
flux cosmological tests (see, e.g., pg. 798]). On the other hand, for advanced detectors 
and interesting cosmological models, the distribution P{p\pQ,C,V) is sensitive to h at the 
10% level, which makes possible a measurement of h from the observations {p|p > po} alone. 
Figure ^ shows P(p|po,X)/P(p|po, Co, 2^) for two Einstein-deSitter cosmological models that 
differ only by h (in both cases Mo = 1.19 M0). The general trend is that with increasing h 
inspiral events are shifted toward larger p compared to the distribution in a static universe. 



2. Examples: Uniformly distributed neutron star masses 

In a less extreme limit assume that neutron star masses are bounded above by and 
below by mi, that between these bounds their masses are uniformly distributed, and that 
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the component masses of an inspiraling binary system are uncorrelated. Then V{Aio) is 
given approximately by equation p. 14 , and 
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(4.24b) 
(4.24c) 

(4.24d) 

(4.24e) 
(4.24f) 
(4.24g) 

(4.24h) 
(4.24i) 
(4.24j) 



Over the range of cxq relevant for neutron star binary inspiral observations in LIGO-like 
interferometric detectors, dN/dt and d^N/dtdp are only weakly dependent on m/ and m/^ 
for constant mo: taking 



(4.25a) 



we find 
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Observations suggest that e is small ||2^, in which case the distribution of p in LIGO obser- 
vations is insensitive to the finite neutron star mass range. 
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F. Chirp mass spectrum P{M.\pQ,T) 



In a homogeneous and isotropic cosmology the rate at which binary inspiral signals 
corresponding to chirp mass Ai and p greater than po ai'e observed is 
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dz (i - kr'^Y''^ 1 + -2 1 + -2 
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(4.26) 



The catalog's a priori chirp mass distribution is thus 

d^N/dtdM 



P{M\po,C,V) 



dN/dt 



(4.27) 



where dN/dt is given by equation [4.16d| and C represents assumptions regarding the cosmol- 
ogy {Hq, qq, Qq), evolution and the distribution of neutron star masses {S, V), and other, 
unenumerated, model assumptions. 

As discussed in section [11 C| , observational and theoretical evidence suggest that the ns-ns 
binary intrinsic chirp mass distribution is sharply peaked. Neglecting evolution and taking 
the intrinsic chirp mass of all binary systems to be a constant Aio, 
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(4.28) 



(4.29) 



Figure ^ shows the distribution P{Ai\po,C,V) for catalogs with p greater than 8 com- 
piled by an advanced LIGO detector (ro = 355 Mpc) in several matter dominated FRW 
cosmological models. The intrinsic chirp mass of all systems is assumed to be 1.19 M0. Six 
different models are shown, exploring two different h (0.5 and 0.8) and three different go 
(1/4, 1/2 and 3/4). The closely spaced curves with co-located extrema are of the same h 
and differ only in go- Note the strong dependence of P{M.\po,C,V) on h and the weaker, 
but still significant, dependence on go: the dotted and solid curves correspond to flat cos- 
mological models (go = 1/2), the long-dash and dot-long-dash curves correspond to open 
cosmological models (go = 1/4), and the short-dash and dot-short-dash curves correspond to 
the closed models (go = 3/4). In general the smaller go, the more compressed the spectrum 
and the smaller the tail at large Ai. 

Since the abscissa is related to redshift according to = (1 + z)A4o figure H also 
shows the redshift of the preponderance of catalog events. For h = 0.8 most events 
redshift of 9%, while for h = 0.5 most events are at a redshift of 14%. 

More generally, neutron star masses are not all identical; correspondingly, V{A4o) is 
not as simple as a 5-function. Modeling the neutron star mass distribution as uniform be- 
tween lower bound m; and upper bound rrih leads to the intrinsic chirp mass distribution 
P(A^o) given by equation p.ll| ; the corresponding spectrum P{A4\pq,C,V) can be deter- 
mined through equations |4.26| and |4.27| . Figure |^ shows the spectrum P{M\pq, C, V) for four 
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different matter-dominated Einstein-deSitter cosmological models corresponding to two dif- 
ferent h (0.5 and 0.8) and two different neutron star mass distributions: one set of curves 
corresponds to the assumption that M.q = 1.19 for all binary systems, while in the 
second set the binary component masses are assumed to be uniformly distributed between 
1.29 Mq and lAbModot. For all models shown tq = 355 Mpc and po = 8. 

As mi and nih approach mo, 'P(A^o) approaches S{A4o — mo). The dependence of 
P{J^\po,I) on qo, shown in figure ^, is similar but not identical to its dependence on nih — mi 
(for constant mo): variations in go shift the spectrum's large Ai tail, while increasing rrih—mi 
increases the spectrum's overall breadth. 

V. CONCLUSIONS 

Observations of binary inspiral in a single interferometric gravitational wave detector 
can be cataloged according to signal strength (as measured by signal-to-noise ratio p) and 
chirp mass A4. The distribution of events in a catalog composed of observations with p 
greater than a threshold po depends on the Hubble expansion, deceleration parameter, and 
cosmological constant, as well as the distribution of component masses in binary systems and 
evolutionary effects (though for neutron star binary observations evolution is not expected 
to be important). In this paper I find general expressions for the distribution with p and 
Ai of cataloged events, valid in any homogeneous and isotropic cosmological model; I also 
evaluate those distributions explicitly for matter-dominated Friedmann-Robertson- Walker 
models and simple models of the neutron star mass distribution. 

These distribution have two immediate, practical uses in gravitational- wave data analysis 
for interferometric detectors: first, when evaluated for with the cosmological parameters 
reflecting our current best understanding of the universe, they are the prior probabilities 
which, together with the matched-filtered detector output, form the likelihood function 
and determine the posterior probability that an inspiral has been detected; second, when 
compared with the observed distribution in p and Ai of many separate binary inspiral 
observations, they are used to infer new and more informed estimates for the cosmological 
parameters that describe the universe. 

The signal-to-noise ratio p of a binary inspiral event is a detector-dependent measure 
of the signal strength from a radiation source that is estimated in the course of making 
an observation. The normalization of p depends on the detector's noise power spectral 
density Sn{f)- In interferometric gravitational wave detectors like LIGO and VIRGO the 
normalization involves a characteristic distance tq and a function of the detector bandwidth 
C, both of which are depend on Sn{f)- 

The characteristic distance tq gives an overall sense of the depth to which the detector can 
"see" binary systems whose radiation traverses the detector bandwidth, which is determined 
by (. Advanced LIGO interferometers are expected to be most sensitive to binary inspiral 
radiation in the bandwidth 30-200 Hz: over 90% of the signal-to-noise ratio is contributed 
by signal power in this narrow band. 

When searching the output of a detector for the signal from a source an accurate model 
of the detector response is needed. The model need not be accurate over the entire frequency 
domain, but only over that part of the domain where the signal power overlaps with the 
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detector bandwidth. The dependence of p on ^ suggests the definition of a bandwidth func- 
tion B for binary inspiral observations that will be useful for determining over what range 
approximate templates describing the detector response need be accurate. 

Cosmological tests based on catalogs of binary inspiral observations with p greater than 
a threshold po depend on the distribution of cataloged events with p and M.. Tests can 
make use of all the information available in the catalog or properly constructed summaries 
of the cataloged events. The sensitivity of the test depends on how the expected catalog 
distributions change with changing cosmological models. 

For advanced LIGO detectors, the most distant neutron star binary inspiral events with 
signal-to-noise ratio greater than 8 will arise from distances not exceeding approximately 
2 Gpc, corresponding to a redshift of 0.48 (0.26) for h = 0.8 (0.5). The depth is only 
weakly dependent on the range of neutron star masses or the deceleration parameter. As 
the binary system mass increases so does the distance it can be seen, up to a limit: in 
a matter dominated Einstein-deSitter cosmological model with h — 0.8 (0.5) the hmit is 
approximately 2; = 2.7 (1.7) for binaries consisting of approximately 10 Mq black holes. 

The distribution of catalog events with p depends primarily on h and is only very weakly 
sensitive to either go or the range of neutron star masses; however, the chirp mass spectrum 
{i.e., the distribution of catalog events with J\4) is very sensitive to h and reasonably sensitive 
to both qq and the neutron star mass range. This suggests that the spectrum is an especially 
powerful tool for cosmological measurements. The dependence of the spectrum on qq and 
the neutron star mass distribution is similar (though not identical), suggesting that it may 
be difficult to determine these separately from observations. I am currently investigating 
this point and will report on it at a later time. 
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FIG. 1. The signal-to- noise ratio of the radiation from an inspirahng compact binary in an 



interferometric gravitational wave depends, through the function ^ (defined in eq. 3.6e ) on the 
redshifted orbital frequency /max of the system's last orbit before coalescence. Here is shown 
C(/max) for initial LIGO (dashed), VIRGO (dotted) and advanced LIGO (solid) interferometers. 
For more detail, see the discussion in section |IIIF| and HID , 
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FIG. 2. The bandwidth function B{f) (defined in eq. ^.lOf ) describes the fraction of con- 
tributed by signal power at frequency / in a logarithmic frequency interval. Here is shown B for 
initial LIGO, VIRGO and advanced LIGO interferometers. The bandwidth function is especially 
useful for determining over what frequency interval an approximate model of the detector response 
to the radiation, which might be used for identifying the presence of a signal in detector output, 
must accurately mimic the real detector response. For more details see the discussion in section 

HID 21. 
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FIG. 3. The distance to the farthest mspirahng binary system with signal-to-noise ratio p 
greater than a threshold po depends on the detector noise spectrum, the binary system component 
masses, and the cosmological model. Shown here is the redshift to the farthest ns-ns binary system 
observable with /) > 8 in an advanced LIGO detector as a function of the Hubble parameter h 
(the Hubble constant in units of lOOKm/s/Mpc). The three curves represent matter dominated 
Friedmann-Robertson- Walker cosmological models with different qo- a closed model with qq = 3/4 
(dashed curve), a flat model {qo = 1/2, solid curve), and an open model with qo = 1/4 (dotted 
curve). For more discussion see gVD| . 
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FIG. 4. The same as figure ^, except shown here is the luminosity distance instead of the 
redshift. 
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FIG. 5. The expected redshift to the farthest inspiraUng binary system observed by an advanced 
LIGO detector with p > 8 in matter-dominated Einstein-deSitter cosmological models (go = 1/2)- 
Results for symmetric binary systems, consisting of two components each of mass Mb^, and asym- 
metric binary systems, consisting of a 1.4 component and a component, are shown for 
h = 0.8 and h = 0.5. The maximum observable depth at any signal-to-noise threshold is limited 
by the cosmological model and the properties of the detector; for advanced LIGO interferometers 
and po = 8 it peaks for symmetric binaries composed of 10 Mq black holes. For more discussion 



see section IV D . 
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FIG. 6. The rate of ns-ns binary inspiral observations with signal-to-noise ratio greater than 
8 in an advanced LIGO detector is largely insensitive to the neutron star mass range or the 
deceleration parameter in matter dominated Friedmann-Robertson- Walker cosmological models. 
The solid curve shows the expected rate in an Einstein-deSitter model as a function of the Hubble 
parameter h assuming the co-moving ns-ns binary coalescence rate density at the current epoch 
is l.l/iMpc~^ yr^^ (solid curve); the dashed curve shows the same assuming the rate density is 
8 X 10~^ Mpc"'^ yr^^, which is independent of h. For more discussion see section IV E| . 
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FIG. 7. The expected distribution of ns-ns inspiral events with p greater than 8 in advanced 
LIGO detectors depends almost exclusively on the Hubble parameter h. Shown here is the ratio of 
the distribution in two matter dominated Friedmann-Robertson- Walker cosmological models to the 



distribution expected in a flat and static cosmological model. For more details see section IV E 1 
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FIG. 8. A binary system's observed chirp mass Ai depends on its redshift; consequently, a 
ns-ns binary inspiral sample will show a range of chirp masses corresponding to the range of 
system redshifts. Shown here is the expected distribution of A4 for binary systems consisting 
of two 1.37 Mq neutron stars with p > 8 m advanced LIGO detectors for open {qq = 1/4), flat 
(go = 1/2) and closed (go = 3/4) matter-dominated Priedmann- Robertson- Walker cosmological 
models with h = 0.5 and h = 0.8. In all cases, as qo increases the tail of the chirp mass spectrum 



is extended. For more details, see the discussion in IV F. 
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FIG. 9. Neutron stars do not all share the same mass, although indications are that the mass 
range is small. Shown here is the expected chirp mass distribution for observations with p > 8 made 
in an advanced LIGO detector for two different matter-dominated Einstein-deSitter cosmological 
models (go = 1/2 and h = 0.5 or h = 0.8) and two different neutron star mass distributions. In 
the first distribution all binaries are assumed to have intrinsic chirp mass Mq = 1.19 Mq, while in 
the second the binary component masses are assumed to be uniformly distributed between lower 
bound 1.29 Mq and upper bound 1.45 Mq. As the mass distribution broadens, the chirp mass 
spectrum also broadens. It does so nearly symmetrically; in contrast, variations in qq for fixed 
mass distribution (shown in figure ^ alter the large M tail of the spectrum, leaving the small Ai 



tail essentially unchanged. For more details see section IVF 
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